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Abstract 

A theory of nonmagnetic impurities in an anisotropic superconductor in- 
cluding the effect of anisotropic (momentum-dependent) impurity scattering 
is given. It is shown that for a strongly anisotropic scattering the reduction 
of the pair-breaking effect of the impurities is large. For a significant over- 
lap between the anisotropy functions of the scattering potential and that of 
the pair potential and for a large amount of anisotropic scattering rate in 
impurity potential the superconductivity becomes robust vis a vis impurity 
concentration. The implications of our result for YBCO high-temperature 
superconductor are discussed. The experimental data of electron irradiation- 
induced Tc suppression (Phys. Rev. B50, 15967 (1994)) is understood quan- 
titatively and a good qualitative agreement with the ion (A^e+) damage and 
Pr substitution-induced Tc decrease data (Phys. Rev. B50, 3266 (1994)) is 
obtained. 
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There now exists a considerable experimental evidence supporting the d-wave supercon- 
ductivity in the cuprates (for review see Refs. 1-4). Nevertheless, this scenario still faces 
some theoretical difficulties. One of these is the predicted extreme suppression of the crit- 
ical temperature Tc by nonmagnetic impurities. Experimentally, however, the observed 
suppression of by impurities or radiation damage in YBCO is much more gradual.01ii 
This issue was critically examined by Radtke et al.i who considered isotropic impurity scat- 
tering within the second Born approximation by applying the Eliashberg formalism. Their 
predictions in both weak- and strong-coupling theory gave a Tc suppression which was close 
to the Abrikosov-Gorkov scaling functionUli with an effective impurity scattering rate. This 
led to an approximate universal dependence of on the planar residual resistivity po, 
which did not depend on the details of the microscopic pairing. In order to verify the re- 
sults of Radtke et al.i systematic electron irradiation experiments on YBCO were carried 
out by Giapintzakis et ahlll The measured initial slope of impurity induced Tc suppres- 
sion was dTc/dp ~ — 0.30 K / pilcm^ whereas the predicted value was in the range from 
~ —0.7AK/ pricm to —1.2i^//if2cm.i While discussing the experimental results in Ref. 12 
the authors invoked the issue of the anisotropic impurity scattering. They understood their 
data within a model of Millis et al.i assuming a value of 0.5 for a dimensionless parameter 
gi which describes the anisotropy of the scattering potential and modifies the bare isotropic 
impurity scattering rate 1/r according to l/r* = (1 — gi) /r, where 1/r* is the effective 
scattering rate. Thus the analysis by Giapintzakis et ahlll brings out the significant role of 
the anisotropic scattering in understanding the impurity effect on rf-wave superconductivity 
and calls for more detailed theoretical studies. 



In this paper we consider in detail the problem of nonmagnetic impurities in an 
anisotropic superconductor for the case of anisotropic (momentum-dependent) impurity scat- 
tering by applying weak-coupling approximation. We find a remarkable change in the Tc 
suppression which becomes more gradual when the anisotropy function defining anisotropy 
of the impurity potential overlaps with the anisotropy function of the order parameter. Al- 
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though our formahsm is general and vahd for any superconducting order parameter described 
by a one dimensional (ID) irreducible representation of the crystal point group we discuss 
the results for a d-wave superconductor in the context of high temperature superconductiv- 
ity. In a certain limit, the effective scattering rate in our model is identical to that of Millis 
et al.i We compute Tc as a function of planar residual resistivity. Within a certain range of 
scattering potential parameters values we find a quantitative agreement of our results with 
the electron irradiation data.lii Also for an appropriate choice of the impurity potential coef- 
ficients a good qualitative fit to the Pr substitution and A^e"^-irradiation dataS is obtained. 
We take h = = I throughout the paper. 



We consider randomly distributed nonmagnetic impurities in an anisotropic supercon- 
ductor. Treating the electron-impurity scattering within second Born approximation and ne- 
glecting the impurity-impurity interaction,Ei the normal and anomalous temperature Green's 
functions averaged over the impurity positions read 

G(..,k) = %±^. (1) 



A(k) 

C^2+^,2^|A(k) 



F {uj, k) = =^^^4 (2) 

^ ' ^ ~9 , ^ 2 , I A /I \ 12 W 



where the renormalized Matsubara frequency Cj (k) and the renormalized order parameter 
A (k) are given by 

(k) = + lUi \w{k- k') {u, k') (3) 
J (2vr) 

A (k) = A (k) +nj\w{k- k') \'F {uj, k') ^ (4) 
J (27r) 

In above uj = 7rT(2n+l) (T is temperature and n is an integer number), is the quasiparticle 

energy, is impurity (defect) concentration, w {k — k') is momentum dependent impurity 

potential and A (k) is the orbital part of a sing let0 superconducting order parameter defined 

as 



A(k) = Ae(k) (5) 

where e (k) is a real basis function of a ID irreducible representation of an appropriate point 
group, which seems to be good approximation for high Tc superconductors.0 We normalize 
e(k) by taking (e^) = 1, where < ... >= Jp^ dSkn (k) (...) denotes the average value over the 
Fermi surface (FS), n (k) is the angle resolved FS density of states normalized to unity, i.e. 
JpgdSkU (k) = 1, and JpgdSk stands for the integration over the Fermi surface. 



The impurity scattering potential is assumed to be separable and given by 

|«;(k-k')P = KP + ki|V(k)/(k') (6) 

where \wq\ {\wi\) is isotropic (anisotropic) scattering amplitude and / (k) is the momentum- 
dependent anisotropy function. We assume that the overall scattering rate is determined by 
the isotropic component and impose the constraints 

kip<kor, (/) = o, (f) = i (7) 

Therefore the Fermi surface average of the scattering potential is (|u' (k — k') p) = |woP 
and the momentum- dependent part in Eq. (^) represents the deviations from the isotropic 
scattering. It is clear that this kind of anisotropic scattering cannot affect the properties of 
the isotropic superconductor, but it can play a certain role in the case of a superconduc- 
tor with an anisotropic order parameter. Although the structure of scattering potential is 
postulated in Eq.(P) this approach is rather general since no additional assumption about 
/ (k) is made in contrast to previous methods.@lll We note from Eq.(^) and from the form of 
impurity potential (Eq. (|^)) that u is k-dependent. This means that the electron self-energy 
due to impurity scattering and consequently the quasiparticle life-time are anisotropic and 
change over the Fermi surface. Further, it yields from Eqs. (^) and (^ that the impurity 
scattering may change the symmetry of the renormalized order parameter A (k) depending 
on the / (k) symmetry. In this respect our approximation differs from that by Markowitz 
and KadanofiEl who assumed only a change of a degree of order parameter anisotropy but 
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not the anisotropy function itself. Moreover in Ref. 15, the anisotropy of the order pa- 
rameter was introduced in a way appropriate for weak anisotropy only. In the more recent 
study of anisotropic scattering by Millis et al.i the authors also assumed that the anisotropic 
impurity potential does not change the symmetry of the electron anomalous self-energy. We 
may also mention that our approach is different than that by Brink and Zuckermann,lli 
where the scattering potential was essentially isotropic but its amplitude varied with the 
superconducting channels. 



To proceed further, we restrict the wave vectors of the electron self-energy and pairing 
potential to the Fermi surface and replace / d^k/ (27r)^ by Nq JpgdSkn(k) J d^k, where Nq 
is the overall density of states at the Fermi surface. Using Eqs. (|l|), (0), (^) and (|) in Eqs. 
(0) and (^) and performing the integration over C,k (particle-hole symmetry of quasiparticle 
spectrum is assumed) we write 



u (k) = u![l + u{uj, k)] 
A (k) = A [e (k) + e (lu, k)] 



(8) 
(9) 



where u (u;, k) and e (w, k) separate into the isotropic (subscript s) and anisotropic (subscript 
a) parts as follows 



u {u, k) = Us (uj) + Ua (uj) f (k) 
e (uj, k) = Cs (uj) + Ca (u;) / (k) 

which are determined by the self-consistent equations 

1 + uiu), k) 



(10) 
(11) 



Us iuj) = To / dSkTi (k) 

JFS 



CU2+ |A(k) |2 



1/2 



u,{uj) = tJ d5'fcn(k)/(k) 

JFS 



1 + u (cu, k) 
'Cj^ + |A (k) |2' 



1/2 



Cs {uj) = To / dSkU (k) — ^ 

•IFS 



e (k) + e (cu, k) 



Cj^ + |A (k) 



1/2 



(12) 



(13) 



(14) 
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e„M =ri / dSkn{k)fik) 

JFS 



e(k) + e(cu,k) 

1/2 



(15) 



CD2+ |A(k) |2 

In writing above we have introduced the isotropic Fq and anisotropic Fi impurity scattering 
rates (Fi < Fq) 



To = 7rNoni\wo\'^, Fi = vrA^o'^ikiP 
The gap function is given by the weak-couphng self-consistent equation 

A (k) = -T V V (k, k') 

with the phenomenological separable pair potential V (k, k') taken as 



(16) 



(17) 



F(k,k') = -^e(k)e(k') (18) 
Following standard procedure,^ we obtain the equation for the critical temperature Tc as 



(19) 



with 



(/(^))a=o = / dSkn(k)4 

JFS OJ 



e(k) 



ujo (k) 



A(k) 



(20) 



A=0 



where is the critical temperature in the absence of impurities and uj^ (k) 
Using Eqs.(|8|)-([T5|), we get for a A ^ hmit 

(Do (k) = + V-sign [oj] 

and 



(D(k) 



A=0- 



(21) 



A(k) 
A 



e(k) + ^(e) + 



\LO\ 



ief) f (k) 



(22) 



jA=o 1^1 l^l+Fo-Fi 

The last two terms on the righthand side of Eq. (^) represent the possible impurity induced 
anisotropy of the renormalized order parameter, which is absent in the works of Markowitz 
and KadanoflJli as well as of Millis et al.i, where A (k) = Ae (k) is assumed. It may 



be mentioned that A (k) and A (k) have the same anisotropy given by e (k) function if 
/ (k) = ±e (k) ((e) = 0) only. Based on Eqs. (|0D, (|l]), and (||) we get from Eq. 



In 



Co 



{{ey 



27lTr 



/ /\2 ^1 



E 

c n>0 



1 



(23) 



n 



2 27TTr 



n 



2 27rT, 



where if) {z) is digamma function.lla The first term on the righthand side of Eq. ( p3D gives 
the Tc suppression due to the isotropic scattering. Since the second term which couples the 
anisotropy functions e (k) and / (k) is always nonnegative, T^. does not decrease as fast as 
for the isotropic scattering only. In other words, an anisotropic potential of the form given 
by Eq. (|^) diminishes the suppression of superconductivity if the scalar product (e/) value 
is nonzero, which may be the case in many cuprate superconducting compounds. For an 
isotropic superconductor e (k) = 1 ((e) = 1, (e/) = 0) and it yields from Eq. (|23|) that the 
critical temperature does not depend on the impurity scattering which is in accordance with 
the Anderson's theorem.lli Finally, Eq. (|2^) may be written in a more compact form as 



CO 



2nTr 



+ 



(24) 



(efY 



27rT, 



Our model has two more dimesionless parameters than the isotropic scattering model. 
First is (e/)^ = [Jpg dSkU (k) e (k) / (k)]^, which describes the interplay between the pair po- 
tential V (k, k') (Eq. (p^)) and the anisotropic part of the scattering potential jti? (k — k') p 
(Eq. (^)). This parameter is determined by the symmetry of the superconducting state 
(e (k)) and that of the impurity scattering matrix element (/ (k)). According to the normal- 
ization of the order parameter orbital function (e^) = 1 (Eq. (^) and that of the anisotropy 
function of the impurity potential (/^) = 1 (Eq. (^) the parameter (e/)^ takes values 
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between and 1. When (e/)^ = then the e (k) and / (k) functions are orthogonal, which 
means that the pair potential V (k, k') and the impurity scattering potential \w (k — k') p do 
not couple and the Tc decrease in Eq. p^) is due to isotropic scattering only. On the other 
hand for (e/)^ = 1 we deal with / (k) = ±e (k) and the pair-breaking effect is minimized by 
the anisotropic part of the scattering potential which is proportional to the pair potential. 
In this case Eq. (|^) yields (e) = and Eq. (^) becomeJ^l 

with gj = Ti/Tq which leads to the pair-breaking parameter (1 — (7/)ro/ (27rTc). On the 
other hand if we calculate gj coefficient defined in Ref. 8 with the impurity potential from 
Eqs. (§|)-(0) we get gi = (e)^ + (e/)^ (Ti/To), which reduces to our value of gi for (e/)^ = 1 
and (e) = 0. Thus in this case our pair- breaking parameter is identical to the one obtained 
by Millis et al.i The second parameter in our model (Ti/Tq) represents the amount of 
anisotropic scattering rate in impurity potential normalized by the isotropic scattering rate 
(Eq. its value ranges also from to 1. For Ti/Tq = we obtain 



which is the isotropic scattering case@ and yields a considerable critical temperature sup- 



pression for (e) 7^ 1. When (e) = then Eq. ( P^ ) gives a suppression curve for a d-wave 
superconductors with isotropic scattering and is the weak-coupling version of the form used 
by Radtke et al.i In the case of strong anisotropic scattering Ti/Tq = 1 and the Tc equation 
reads 

It is easy to see that the critical temperature suppression becomes more gradual now and 
may be even reversed into Tc increase for a significant overlap between e (k) and / (k) func- 
tions, that is when (e/)^ ~ 1. 
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Our results for the dependence of Tc/Tc^ on the isotropic scattering rate Tq/2hTcq are 
shown in Figs. la-Id for a selection of the model parameters (e/)^ = 0.2, 0.4, 0.8 and 0.95 
and Fi/Fo = 0, 0.5, 0.9, 0.95 and 1.0. We have assumed here (e) = 0.0 Based on these we 
make the following remarks: (1) In all curves the depression of Tc in the limit of impurity 
concentration — > is given by the initial slope 



Fn 



J- 



(28) 



which decreases drastically as (e/)^ Fi/Fq approaches unity; (2) For a given value of Fi/Fq, 
the value of Tq/2'kTc^ needed to suppress superconductivity increases as (e/)^ is increased; 
(3) When there is a significant overlap between the anisotropy functions e (k) and / (k), e.g. 
(e/)^ ~ 0.8 (Fig. Ic) the value of needed to destroy superconductivity is increased 

considerably when Fi/Fq becomes large. 



In order to make contact with experiment, we estimate the planar residual resistivity 
Po, which is a normal state property and according to Eq. (PT| ) depends on the isotropic 
scattering rate Fq exclusively. It is worth mentionnig here that in the normal state the 
influence of the impurity scattering on the electron self-energy is reflected by the frequency 
rescaling only (Eq. (^)), and hence the scattering process is characterized by Fq parameter 
completely. Therefore neither of anisotropic scattering parameters enters the equations 
determining the normal state properties. Using a Drude form of the low frequency residual 
electrical conductivity at zero frequency a = Upir/AiT, where Upi is the plasma frequency 
and 1/r = 2Fo, we representil the planar residual resistivity in terms of the dimensionless 
pair-breaking parameter Fo/27rTcQ 

Po ^ 10.18 X 10-2 X ^T,„ ( ) /^^ cm (29) 

^pi V^T^oJ 

with ujpi in eV and Tc^ in Kelvin. From Eqs. (^) and (^) we get the initial slope for a 
(i-wave superconductor 
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dpo 



0.615 X col, - (e/)' T^) ^V/"^ 



cm 



(30) 



In the electron irradiation experiment in YBCO Giapintzakis et al.c3 obtained dTc/dp ~ 
—0.30 ± 0.04 K/fiQ cm (p is resistivity at 145 K and dTc/dp ~ dTc/dpo). Taking 
the plasma frequency Upi ranging from 1.1 to 1.4 eV, which is the experimental esti- 



reproduced by the anisotropic scattering parameters with values given by a constraint 



uncertainty of the plasma frequency and the dTc/dp measurement accuracy.EJ Our calcula- 
tion focused entirely on a single Cu02 plane seems to be a good approximation here since 
the low-energy electron irradiation, used in this experiment, displaces the oxygen atoms only 
and an appropriate measurement method probes the contribution to suppression of only 
the oxygen defects on the Cu02 planes.lii The two dimensional approach is not so justified in 
the interpretation of the experimental data of Ref. 11 where Pr substitution and ion (A^e"*") 
damage were applied. The Pr substitutes onto the Y site and a similar defect is probably 
induced by ion irradiation since the suppression induced by both methods is analogous. 
Considering this caveat, our theoretical results and the experimental data of Ref. 11 are 
shown in Fig. 2 for an illustrative purpose mainly. The data were read from Fig. 4 of Sun 
et al.0 and the region between the curves corresponds to the Tc computed from Eqs. ( p^ 
and ( p9D with (e) = 0, (e/)^ = 0.95, and Fi/Fo = 0.96 for plasma frequencies Upi ranging 
from 1.1 to 1.4 eV. We did not try to adjust the amount of anisotropic scattering present 
in our model so as to get a best fit to the data. Nevertheless, we note a good qualitative 
agreement of the theoretical results with the experimental data of Ref. 11. The experimen- 
tal data show, however a long tail Tc suppression which is not refiected in the computed T^. 
We think that this feature may be due to a slight orthorhombic anisotropy of the system,0 
which was neglected in our calculation by the assumption of (e) = 0. 

Before concluding, we give some critical remarks concerning our approach. We have 



mate 




experimental data can be 
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employed a weak-coupling approximation neglecting the strong-coupling corrections. We 
expect that as in Ref. 6, the strong-coupling effects would rescale the scattering rates. Fur- 
ther, we have neglected the interaction between the nearest Cu02 planes, restricting our 
considerations to a single copper-oxide plane. This simplification may not be valid for the 
interpretation of the experimental data of Ref. 11, where the defects are not in the Cu02 
planes. Finally, we have assumed a model separable momentum-dependent impurity po- 
tential, which is obviously not the most general way of treating the problem, but is more 
general than the one applied in the previous studies.i'El 



In summary, we have given a theory of anisotropic impurity scattering in anisotropic 
superconductors. The impurity potential is assumed to be separable according to Eq. (^. 
There are two parameters characterizing the scattering anisotropy in our approach. The 
first of them ((e/)^) represents the interplay between the symmetry of the superconducting 
order parameter and that of the impurity potential, the second (Fi/Fo) gives the amount 
of anisotropic scattering versus the isotropic one. We find that for a significant overlap 
between the pair potential and the impurity potential that is for large (e/)^ values, and 
for a large value of Fi/Fq, the anisotropic superconductivity becomes robust vis a vis the 
impurity concentration. The experimental data of the electron irradiation-induced sup- 
pression m YBColi is understood quantitatively within our model. We also obtain a good 
qualitative agreement with the observed T^. decrease in YBCO due to ion (A^e"*") damage 
and Pr substitution.0 



We thank J. Taylor for a numerical assistance. This work was supported by the Natural 
Sciences and Engineering Research Council of Canada. 
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FIGURE CAPTIONS 



Fig. 1. Normalized critical temperature Tc/Tcq as a function of the normalized isotropic scat- 
tering rate ro/27rTcQ for different values of the normalized anisotropic scattering rate Ti/Vq = 
0.5 (dotted curve), 0.9 (short-dashed curve), 0.95 (long-dashed curve), 1.0 (dot-dashed 
curve). The solid curve represents the isotropic scattering pair- breaking effect (Fi/Fo = 0). 
We have taken {eff = 0.2 (a), 0.4 (b), 0.8 (c), 0.95 (d), and (e) = 0. 

Fig. 2. Critical temperature Tc (area between solid curves) of a superconductor with the 
critical temperature in the absence of impurities Teg = 90 X vs residual resistivity for 
(e/)^ = 0.95, Fi/Fo = 0.96, (e) = and plasma frequencies cupi between 1.1 and 1.4 eV. 
The experimental data of Ref. 11 for YBCO are shown with circles (Pr substitution) and 
crosses (ion damage). 
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